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' Abstract. An integral relation is established between the Green functions corresponding to two Hamiltonians 

which are supersymmetric (SUSY) partners and in general may possess both discrete and continuous spectra. 
It is shown that when the continuous spectrum is present the trace of the difference of the Green functions for 
SUSY partners is a finite quantity which may or may not be equal to zero despite the divergence of the traces 
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' of each Green function. Our findings are illustrated by using the free particle example considered both on the 



whole real line and on a half line. 
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1. Introduction 



cr. 

At present there is a growing interest in the study of different properties transformations induced 
by supersymmetry (SUSY) in Quantum Mechanics. Recently a special issue of Journal of 
Physics A (see vol. 37, No 43, 2004) was devoted to research work in this subject. Despite of 
the growing number of papers in this field many questions still remain open and require further 
study In particular, the authors are aware of only one paper [T] devoted to the study SUSY 
transformations at the level of Green functions. For the case of a transformation deleting the 
ground state of the initial Hamiltonian, Sukumar has studied an integral relation between the 
Green functions for SUSY partners and has formulated conditions leading to the vanishing of 
some matrix elements of a Hamiltonian and related this property to a hidden supersymmetry 
of the system. Transformation of Green functions is not explicitly discussed in that paper. 
Moreover, we have found that formula (28) of pQ relating integrals over Green functions for 
SUSY partners may need to be corrected if a continuous spectrum is present. 

In this paper we give a simple formula for the Green function of the SUSY partner 
Hamiltonian both for confining and for scattering potentials and generalize results of the paper 
[1J to the case where the continuous spectrum is present. As an application of this general 
formula we consider the case of the Schrodinger equation with a scattering potential defined 
both on the whole real axis and on a half line when the Schrodinger equation is reduced to a 
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singular Sturm-Liouville problem. Regular Sturm-Liouville problem is considered in a separate 
publication [2|. 

2. Green function of the Schrodinger equation 

In this section we cite some properties of the Green function of the one-dimensional Schrodinger 
equation for a spectral problem on the whole real line (see e.g. [HI El) which are useful. 

We consider the Schrodinger equation 

(ho - E)ip = ho = -d 2 /dx 2 + V {x) x e (a, b) (1) 

supplemented by the boundary conditions ip(a) = ip(b) = 0. We will concentrate mostly on two 
cases; these are the whole real line a = — oo and b = oo and the half line a = and b = oo. 
We assume that the spectral set spech of this problem consists of M discrete points with the 
possibilities M = or M = oo and possibly a continuum part filling the positive semiaxis. 

The definition of the Green function used by different authors I21I3JE] ma y differ by a constant 
factor. We use a definition of the Green function represented as the kernel of the operator 
(ho — E)~ l as an operator defined in the corresponding Hilbert space. It is well defined for 
all E ^ spec/io- It has two different but equivalent representations. The first representation is 
obtained with the help of two real solutions of equation ((TJ) with a fixed value of the parameter 
E £ spec/io, fio an d fro ("left" and "right" solutions), satisfying zero boundary conditions: 
fio( a ) — 0, fro(b) = 0. Since they correspond to the same E their Wronskian W = W(f r0 , fi ) 
does not depend on x and is a function of E only, and the Green function is 

G (x, y, E) = f l0 (x, E)f r0 (y, E)/W x<y (2) 

G (y,x,E) = G (x,y,E). (3) 
These formulae are clearly equivalent to 

G (x, y, E) = [f l0 (x, E)f r0 (y, E)Q(y - x) + f l0 (y, E)f r0 (x, E)Q(x - y)]/W (4) 

where is the Heaviside step function. 

If the operator h is essentially self-adjoint the set of its discrete spectrum (if present) 
eigenfunctions {ip n }, n — 0,1, . . . ,M, (VVilVw) = &nm together with the continuous spectrum 
eigenfunctions (also if present) ipk, E = k 2 > 0, (^fclVv) — <^(^ — k'), (ip n \ipk) = is complete 
in the Hilbert space 

M r 

^2^n(x)^* n (y) + / dkip k (x)ipl(y) = 5(x-y) 

n=0 J 

and the second representation of the Green function may be found in terms of this set as follows: 

For the spectral problem on the whole real axis the continuous spectrum is two-fold degenerate 
and the integrals over k run from minus infinity to infinity and for the problem on a half line 
they run from zero to infinity. 
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It is well-known (see e.g. jH]) that there exist three kinds of SUSY transformations: 

(i) deleting the ground state level of h 

(ii) creating a new ground state level and 

(iii) purely isospectral transformation. 

In all cases the partner Hamiltonian hi = —d 2 /dx 2 + V\ for h Q is defined by the potential 

Vi( x ) — Vo( x ) ~ 2w'(x) w(x) = [log (6) 

where u is a real solution to the initial equation (ho — a)u = with a known as the factorization 
constant. We adopt the notation that a derivative with respect to x is denoted by the prime 
symbol. To provide a nonsingular potential difference a should be less than or equal to the 
ground state energy of ho if it has a discrete spectrum or lower than the continuum threshold 
otherwise. The functions <p n = Lip n , n — 1, 2, . . . , M describe (unnormalized) bound states and 
ifE = Lip E correspond to (unnormalized) scattering states of h\. Here 

L = —dj 'dx + w(x) (7) 

is the transformation operator (intertwiner) satisfying Lh = h\L and Lu = 0. The 
normalization constants are easily obtained with the help of the factorization property L + L = 
h Q — a where L + = d/dx + w(x). The functions 

Xn = (E n - a)- 1 ' 2 ^ Xe = (E- a)- l / 2 L^j E (8) 

form an ortho normal set. 

Theorem 1 Let Gq(x, y, E) be the Green function for ho. Then for all three cases enumerated 
above the Green function for h% is 



Gi(x, y, E) = — [L x L y G (x, y,E)-5(x- y)\ . (9) 
hi — ct 



In case (ii) it has a simple pole at E = a. In cases (i) and (iii) it is regular at E = a and can 
be calculated as follows: 



Gi(x,y,a) 



dGp(x ) y 1 E) 
dE 



(10) 



E=a 

Here L x is the operator given in |7|) and L y is the same operator where x is replaced by y. 
Proof. In case (i) u = ipo and the set {xe, Xn, n — 1,2, ... , M} is complete. Therefore 

G^y,E)=±^^ + f^Mdk. (11) 

71=1 11 

Now we replace x using (jHJ) which yields 

Gi(x,y,E) = 

l^L x L y (Ef=ife " E^slM^mv) + fdki^ ~ j^bIMWM)- 
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The statement for E ^ a follows from here if in the first sum and in the first integral we express 
Lip in terms of Xi make use of the completeness condition for the set \ and formula (0) for Gq. 
The fact that here the sum starts from n = 1 and in (j3J) it starts from n = cannot cause any 
problems since Lipn = 0. For E = a formula (JTTjl can be written in the form 



Gt(x,y,a) 



M 



—l l ( v t^^Ml + f M^MM* 

9E x v \J^ E n — E J k 2 -E 



E=a 



from which (|10|) follows in this case. 

In case (ii) let Xa ~ 1/w be the normalized ground state function of h\ corresponding to the 
new discrete level E = a. Then 

Gl (x, y,E) = Y + XMXM + [ X k (x)xt(y) dk {u) 

^ E n - E a- E J k 2 - E 

Now the use of exactly the same transformations as in case (i) reduces f)12|) to Q. 

In case (iii) we start from the same formula with the only difference that the sum now 
starts from n = and following the same line of reasoning as in the earlier cases we get formula 
(jHJ). It is interesting to notice the intermediate result 

Gi(x, y, E) = — l — L x L y [G (x, y, a) - G (x, y, E)\ (13) 
a — hi 

which makes clear how formula (fTUj) arises for this case by taking the limit E — > a. The fact 
that in case (ii) the function Q has a simple pole at E = a is a consequence of the equivalence 
between (JUJ) and ©. □ 

Corollary 1 In terms of the special solutions fio and f r o of the Schrddinger equation for ho 
the Green function G\ for all three cases listed above may be expressed as follows: 

1 

(£7-~a)Wo' 

+ 0(x - y)L v f m (v, E)L x f r0 (x, E)] . (14) 

In case (ii) this function has a simple pole at E = a. In cases (i) and (iii) it is regular at 
E = a and can be calculated as follows: 

Gx(x,y,a) 

d Q(y - x)L x f m {x, E)L y f r0 (y, E) + Q(x - y)L y f l0 (y, E)L x f r0 (x, E)~ 



Gi(x,y,E) = — — ^^r[Q(y - x)L x f m (x,E)Lyf r0 (y,E) 



dE Wn 



(15) 



E=a 



To prove these formulae we substitute G as given in (jlj) into and (fTUj) . Taking the derivative 
of the theta functions in (J3J gives rise to the Dirac delta function which cancels out the 
delta function present in ©. Formula (fTlj) is clearly valid since the SUSY transformations 
necessarily preserve the boundary conditions for all E except perhaps for E = a. This 
implies that fn = Lfi vanishes at x = a and f r i = Lf r0 vanishes at x = b. The 
denominator in (|14j) is just the Wronskian of f r i and fn and may be given in the form 
W(f rl J ll ) = (E-a)W(f r0 ,fio) = (E-a)W . ' ' □ 
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The trace of the Green function defined as [3] f G(x, x, E)dx is usually divergent if the system 
has a continuous spectrum. It is remarkable that the trace of the difference Gq(x, x, E) — 
Gi(x,x, E) is a finite quantity which may or may not be equal to zero. In some cases this 
fact may be explained by another remarkable property. It may happen that the difference 
of infinite normalizations (they diverge as 8{x — y) when y — > x) of the continuous spectrum 
eigenfunctions of the two SUSY partners is a finite quantity. 

Theorem 2 Let fi (x,E) and f r0 (x,E) be solutions of the Schrddinger equation for h 
satisfying the zero boundary conditions at the left and right bound of the interval (a, b) 
respectively, and 

f a (x,E) = Lf l0 (x,E) f rl (x,E)=Lf r0 {x,E) (16) 

be similar solutions for hi related with h by a SUSY transformation with a being the 
factorization constant. Let W be the Wronskian of f r0 and f\o, W = W(f r0 , // ) ■ Then 

JjG {x,x,E) - G l (x,x,E)]dx = w ^ ] _ - (17) 
where Q{E) can be calculated by one of the following formulae: 

Q{E) = (frofll)x=b — (frofll)x=a = (flofrl)x=b ~ (flofrl)x=a (18) 
= -W + (flofrl)x=b ~ (frofll)x=a = W + (f r ofll)x=b ~ (flofrl)x=a ■ (19) 

Proof. From Corollary 1 it follows that G\(x,x,E) = Wo ^_ a ) 

Lf l0 (x)Lf r0 (x). While 

integrating this expression over the interval (a, b) one can transfer the derivative present in 
L either from fi to f r o or from f rQ to fio which leads to one of the following integrands 
fio(x)L + Lf r o(x) or f r0 (x)L + Lfio(x). In both cases the factorization property may be used to 
reduce the integrand to (E — a) fio(x) f r o(x) . Thus we arrive at the relation 

/ Gx{x,x)dx = / f lo (x)fro(x)dx- ®^ (20) 

J a W J a W {E - a) 

where Q(E) is given by IjlKj) . To prove ()19j) it is sufficient to notice that 

f l0 (x,E)f rl (x,E) - f r0 (x,E)f n (x,E) = W (21) 

which is a consequence of (116(1 . The identification of the integrand on the right hand side of 
(J20|) as WoGo(x,x) then leads to the result given in (fT7(l . □ 

Using the first of equalities (fTH|) one can rewrite (|T7jl as follows: 

f \r~< I 7m r* I 77i M J 1 i {flofrl)x=b {frofll)x=a /oo\ 

/ [Go(x,x,E)-Gx(x,x,E)]dx= - + — — r (22) 

Ja a- E W {E - a) 

Now for the case (i) where a — Eq if we compare this result with the corresponding difference 
which can be obtained directly from the expressions for Go given by (0) and for G\ given by 
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(II lj) the following feature may be noted: the first term on the right hand side of (|22j) arises from 
the contribution to Green functions from the discrete spectra and the second term, which as 
we show below may be different of zero, is due to the presence of the continuous spectra. Just 
this contribution was neglected in pQ. So, Theorem 2 presents a generalization of the result 
obtained in pQ to the case where a continuous spectrum may be present. As an application of 
this theorem we are going to consider two particular cases of scattering potentials defined both 
on the whole real line and on a semiaxis. 

Corollary 2 If h is a scattering Hamiltonian with the potential V satisfying for the spectral 
problem on the whole line the condition 



[1 + \x\)\Vo(x)\dx < oo 

then for E ^ a, lm\[E > the following equality 

[G (x, x, E) - G x (x, x, E)]dx = ^—^ (23) 



CO 



k z + tan k z + a 1 



holds, where E = k 2 , a = —a 2 ; 5=1 for the case (i), 5 = — 1 for the case (ii) and 5 = for 
the case (Hi). 

Proof. The statement readily follows from the fact that any scattering potential has a pair 
of solutions (Jost solutions, see e.g. jjj) with the following asymptotics at the right infinity 

fl, r (x, E) -> e TiKX E = k 2 Im/t > x -> oo 

and similar asymptotics at the left infinity and the use of an appropriate part of equalities (fTSj) 
and (|inj). The Wronskian Wo for Jost solutions can easily be calculated, Wo = -2k. □ 

So, we see that despite the fact that for both h and hi the continuous spectrum 
eigenfunctions are normalized to the Dirac delta function, (i.e.) that in both cases they have 
equal infinite norms, the difference of these infinities is a finite non-zero quantity in cases (i) 
and (ii) and it is zero in case (iii). 

For instance in case (ii) the following equality arises: 



where 



P(k)= / l\Mx)\ 2 ~ \Xk(x)\ 2 }dx . (25) 



Equation (|2"3j) may be reduced to the Stieltjes transform and the function P may be found by 
the Stieltjes inversion formula (see e.g. jlj). To establish this we first notice that the integral on 
the left hand side of (J2"2j) is different from zero only if P(k) is an even function which we assume 
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to be the case. Therefore it can be considered only for positive ks and we can let k 2 = A. So, 
takes the form 

^ = R(E) 

where the measure p(A) is continuous for A > 0, dp(X) = -^P(X)dX and such that for negative 
As the integral is zero. Now the Stieltjes inversion formula yields 

^ = ^L^ mE )-R(E)] E = X + ir 

where the bar over E denotes the complex conjugate to E. Note that because of the condition 
lm\/E > the square root of E has different signs for E in the upper and lower halves of 
the complex E-pl&ne. Therefore the function R(E) has a cut along the real axis and the jump 
across this cut defines the function P(X). After a simple calculation one gets 

P(X) = an-\X 2 + a 2 )' 1 . (26) 

It must be noted that in the present case the interchange of the integrals over the space variable 
x taken in the difference of (0) and (fTTj) at y — x with the integral over the momentum k is 
justified. 

Corollary 3 If h Q is a scattering Hamiltonian with the potential Vq for the spectral problem on 
a half line satisfying the condition 

x\Vo(x)\dx < oo 



o 



then there exist only two kinds of SUSY transformations keeping the zero boundary condition 
at the origin. If ho has the discrete spectrum its ground state ipo may be deleted, (u = ipo, case 
(%)) and there is a possibility to keep the spectrum unchanged (case (Hi)). The last possibility 
may be realized with u(x) = fi(x, E), E < where fi is such that fi(0, E) = 0. In this case the 
following trace formula is valid: 

/ [G (x, x, E) - s, E)}dx = - 2 , - -^- 2 (27) 

where E = k 2 and a = —a 2 ; for the case (i) 8\ — 62 — 1 and for the case (Hi) 5\ = 0, 62 = — 1. 

Proof.. The proof is based on the fact that for such potentials the left solution goes to zero 
like x when x — > and the right solution has the asymptotics f r Q ~ exp(mx), E = k 2 , Imk > 
(see e.g. jH|). It follows from here that when E is not a spectral point the left solution has a 
growing asymptotics at infinity, f i ~ — ^ exp(— inx\ x — > 00, where Wo is the Wronskian of 
f r o and fio. It is not possible to create a new bound state in this case since the Schrodinger 
equation with such a potential has no solutions going to infinity as x approaches the origin. □ 

Example 1. Free motion on the line, Vq(x) = 0, x 6 R 
The Green function is 

G {x, y, E) = -L e iKlx ~ vl lmn>0 E = k 2 . 

2k 
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The choice u = cosh(ax), a = —a 2 leads to the one soliton potential V\ = — 2a 2 sech 2 (aa;) with 
the Green function 

G 1 {x,y,E) = if K (x)f_ K (y)/lMK 2 + a 2 )] x < y 

where f K (x) = exp(— inx)(in + atanhax), which clearly has a pole at the ground state energy 
E = —a 2 . The residue at the pole is ip (x)ip (y) where ipo(x) = sech(ax) which is just 

the ground state of the one-soliton potential. 

Continuous spectrum eigenf unctions of ho, iftk( x ) = l/v^rexp (ikx) are transformed into 
continuous spectrum eigenfunctions for hi, — [—ik + tanh(ax)] exp (ikx)/ ^2ix(k 2 + a 2 ). 

The direct calculation of the function P{k) given in (j2H) gives exactly the result pHjl. 

Example 2. Free motion on a half-line with zero angular momentum, Vq(x) = 0, x G M + . 
The Green function is 

Go(x, y,E) = — sm(Kx) expUny) E = k 2 \mn > x < y . 

K 

It must be noted that the use of the transformation function u = cosh (ax) gives the same 
one-soliton potential as in the previous example which when considered on a half-line has 
only a continuous spectrum but its continuous spectrum eigenfunctions cannot be obtained by 
applying the transformation operator L to the free particle eigenfunctions since it does not 
preserve the zero boundary condition at the origin. Another choice u = sinh(ax) creates the 
potential V\ = 2a 2 csch 2 (ax) which is singular at the origin and has only a continuous spectrum 
and which may be obtained with the help of the SUSY transformation from the free particle 
Hamiltonian. Its Green function is 



Hi 



my 



G\ (x, y, E) = — [k + iacth(ay)] [k cos(kx) — acth(ax) sm(Kx)} 

k(k 2 + a 2 ) 

E = k 2 Imn > . 

It is not difficult to see that despite the presence of the denominator, G\ is regular for all 
k including the point k = ia and discontinuous along the positive real axis in the complex 
i^-plane. The jump across this cut is proportional to the product of continuous spectrum 
eigenfunctions of h\ which are given by 



= a / ^-[— k cos(kx) + acth(ax) sin (fee)] . 

W 7i(k z + a 2 ) 

In contrast to the previous example the difference — Q. now oscillates when x — >• oo, the 
Riemann integral of this difference over the space variable is divergent and the interchange of 
the integrals over k and x is impossible. Nevertheless, if one assumes that the improper integral 
over x is a limit of a proper integral one can write 

r™PpA)dk = R(E)= 1 (28) 

a^oo n k z — E Ziian — k z ) 
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where E = n 2 , a = —a 2 and P(k, A) = / [xpl(x) — Xk( x )]d% ■ 

Jo 

In our example the function P{k } A) is given by 

2acth(aA)sm 2 (kA) - A;sin(2Afc) 
F[K > A) ~ 7r(k 2 + a 2 ) 

which when substituted into the left hand side of (J2~%|) gives exactly the function R(E). So this 
example shows that in contrast to the previous example where the difference of normalisations 
is a finite quantity, here this value is undetermined. Nevertheless, the contribution to the trace 
of the difference Gq — G\ from the continuous spectra of ho and hi is well defined. 



5. Conclusions 



In this paper we have studied the relation between the Green functions corresponding to two 
Hamiltonians which are SUSY partners. We have shown that it is possible to establish a 
relation between the traces of the Green functions for the two partner Hamiltonians for the 
cases of deletion of the ground state, the addition of a new ground state and when the two 
Hamiltonians are isospectral. The formulae derived in this paper are valid for the general case 
of Hamiltonians having both discrete and continuous spectra. Our results show that when a 
continuous spectrum is present, each of the traces of the Green functions for the SUSY partners 
may diverge but the difference between the traces can be finite. We have illustrated our results 
by considering the case of the free motion on the full line and the case of the free motion of a 
particle with zero angular momentum on the half-line. 

Finally we would like to note that the difference of the traces of the Green functions of 
the two SUSY partner Hamiltonians appears as the trace (actually super-trace) of the Green 
function of the supersymmetric Schrodinger equation (supersymmetric Green function). Thus, 
our results reveal the possibility of divergence of the component traces of the supersymmetric 
Green function while its super-trace remains finite. 
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